Non-reciprocal light diffraction by a vortex magnetic particle 
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We report a theoretical study of light diffraction by a spherical magnetic particle with a vortex 
magnetization distribution. It is shown that the intensity of the diffracted light involves a non- 
reciprocal contribution. This contribution depends on the vorticity of particle magnetization. It 
appears due to the excitation of an electric quadrupole, magnetic dipole and the addition to the 
electric dipole moment in the particle, that depend on the particle magnetization vorticity. The 
estimation of the non-reciprocal contrbution for a cobalt particle and two linear polarizations of the 
incident light fits the data of recent experimental studies in the lattice of triangle magnetic particles 
by an order of magnitude. 

PACS numbers: 75.50.Cc, 75.75.Fk, 78.67.Bf 



I. INTRODUCTION 

Remarkable achievements in micro- and nano- tech- 
nologies open new possibilities in fabrication of opti- 
cal media with novel interesting properties. Photonic 
three-dimensional crystals exhibit extraordinary optical 
effects, such as photonic band gap^ 2 -, negative refrac- 
tive index^, and ultraslow light propagation 5 . Sur- 
face plasmons in planar metal structures lead to en- 
hanced magneto-optical effects 6 - - — , extraordinary light 
transmission through subwavelength holes^, and effec- 
tive generation of the second harmoni o 10 ' 11 . Interest- 
ing objects actively investigated in the past decade are 
planar chiral structures^ - —. Such a structure is ordi- 
narily a two-dimensional regular lattice of non-magnetic 
particles that possesses the time and spatial inversion 
symmetries^!. Nevertheless, light transmission through 
a planar chiral structure experiences polarization rota- 
tion which is different for the light traveling in opposite 
directions due to the absence of mirror symmetry and 
strong anisotropy in the orientation of the structure ele- 
ments. The term "non-reciprocal" is often used for this 
effect, and some speculations appear on the time-reversal 
symmetry violation^ - — . Note that the intensity of light 
diffracted by a planar chiral structure does not change 
with a reversal of the propagation direction, in accor- 
dance with the reciprocity law. The non-reciprocal in- 
tensity effects can appear only in the presence of a mag- 
netic fielder— or in magnetic systems^ - — . Such effects 
have been observed recently in a two-dimensional lattice 
of magnetic vortices 25 . The particles had a triangular 
shape, which allowed manipulation of their vorticity by 
an uniform external magnetic field, thus making it iden- 
tical for all particles. 

The problem of scattering of electromagnetic waves 
by a particle of an arbitrary size with gyrotropic dielec- 
tric permittivity and magnetic permeability has been ad- 
dressed in a general form 2 ! and for magnetic particles 
specifically 2 -' 29 . However, the gyrotropy terms are al- 
ways supposed to be constant over the particle, which 
corresponds to uniform magnetization. This paper is de- 
voted to theoretical calculations of the diffraction of light 



by a magnetic vortex that turns out to be a non-uniform 
magnetic particle. Our theory describes the mechanisms 
underlying the effect of non-reciprocal light scattering 
by such a particle and also allows us to estimate its 
magnitude and compare the theoretical data with the 
experiment^.. Note that this effect is assumed to be sim- 
ply summed over the lattice of particles (i.e., the collec- 
tive phenomena are neglected), thus making it possible 
to focus on the light diffraction by a single particle. 

In Section|TT]we consider some phenomenological argu- 
ments in favor of the non-reciprocal light diffraction by a 
vortex particle. Then a microscopic model of the effect is 
proposed in the approximation of a spherical particle that 
is small compared to the wavelength. The main assump- 
tions are outlined in Section Um Further, we consider two 
approaches to the problem. A simple Born approxima- 
tion (permittivity is close to unity) that helps reveal the 
existence of the effect is described in Section IIV1 Yet, 
within this approximation the effect occurs only for one 
linear polarization of the incident light. Section [V] is de- 
voted to calculation for the arbitrary permittivity, based 
on the first-order perturbation theory with respect to the 
particle size-to- wavelength ratio. We show that the non- 
reciprocal term appears due to excitation of the electric 
quadrupole, magnetic dipole and a small addition to the 
electric dipole moment, that emit light interfering with 
the main magnetization-independent electro-dipole radi- 
aton. Finally, we perform simple estimations for the pa- 
rameters of cobalt and an appropriate particle size and 
compare our calculations to the experimental results 26 . 
A summary of our results is given in Section IVII 



II. PHENOMENOLOGICAL CONSIDERATIONS 

We begin with some phenomenological arguments in 
favor of the non-reciprocal light diffraction by the vor- 
tex particle^ 6 .. If we consider the light scattering cross- 
section summed over the polarizations of incident and 
diffracted light, the reciprocity law takes a simple form 



a (k, k', M (r)) = a (-k', -k, -M (r)) 



(1) 
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here a is the differential cross-section for the diffracted 
light, k and k' are the wave vectors of the incident and 
diffracted beams, M (r) is the magnetization spatial dis- 
tribution. The term "non-reciprocal effect" implies the 
following inequality: 

a(k,k',M(r))^ ( 7(-k',-k,M(r)), (2) 

which can be transformed using Eq. (fTJ) to 

f r(k,k',M(r))^a(k,k',-M(r)), (3) 

thus making it possible to observe the effect by invert- 
ing the magnetization direction instead of swapping the 
source and the detector—. 

For systems without center of inversion the scattering 
cross-section may contain the term ((k + k') • C), where 
C is vector. Being linear in the wavevector, this term 
leads to the non-reciprocal effects described by ((21 EJ • 
The C vector should be a polar vector that also changes 
its sign under the time reversal. For a magnetic scat- 
terer of a centrosymmetrical shape and material C can 
be chosen in the simplest form C = a ([r x M (r)]) which 
is a toroidal moment of the particle associated with the 
magnetic vorticity^, the brackets mean the spatial aver- 
aging over the scatterer, a is a constant. According to 
these considerations, the diffraction of unpolarized light 
by a particle with the vortex magnetization distribution 
is non-reciprocal and the scattering cross-section has a 
contribution depending on the vorticity: 

a (k, k', M (r)) = ... + a ((k + k') • ([r x M (r)]» . (4) 

The weak point of this phenomenological consideration 
is the polarization dependence of the non-reciprocal ef- 
fect. Indeed, the existence of the effect for unpolarized 
light means it exists for at least one linear polarization, 
but the question about the contribution of different polar- 
izations to it arises. However, the experiments currently 
are performed exactly for a linearly polarized lights. 

III. MAIN ASSUMPTIONS AND DEFINITIONS 

As has already been mentioned, we assume that the 
diffracting particle has a spherical shape with radius a. 
The incident wave is assumed to be monochromatic (all 
fields change in time as e 1 "', uj is the wave circular fre- 
quency) and plane. Its wavelength A is much bigger than 
a: 

A >> a. (5) 

The particle has a vortex magnetic moment (see Fig. [T]) 

M = = e y cos0 — e x sin</>. (6) 

Here is the unit vector in spherical coordinates, e x and 
e y are the unit vectors in Cartesian coordinates, <f> is a 
spherical coordinate, magnetic moment M is normalized 
to unity. 




(a) (b) 

FIG. 1. (Color online) The geometry of light diffraction by a 
vortex magnetic particle, (a) The incident light is s-polarized. 
(b) The incident light is p-polarized. 

The magnetic permeability tensor is unit, while the 
dielectric permittivity tensor is assumed to have a locally 
gyrotropic form: 

(e) j7 = sSji + ijeji m M m , (7) 

where Sji is the Kronecker delta, eji m is the completely 
antisymmetric tensor (the Levi-Civita tensor). Since M 
is normalized to unity, the 7 coefficient includes its mag- 
nitude. Expression (JT)) is valid if the magnitude of quasi- 
classic electron oscillations is small compared to the scale 
of the magnetic moment variation. In our case it is equal 
to the particle radius a, hence, we get applicability crite- 
rion 



where e is the absolute electron charge, m is the electron 
mass, Eq is the magnitude of the electric field of the 
incident wave, c stands for the light velocity. The first 
term in the left-hand part of (JSJ is the ratio of the energy 
the electron gains while oscillating in the electric field of 
the wave, to its rest energy. In the optical range this 
ratio is very small for the existing sources, which allows 
simultaneous fulfillment of conditions (|5|) and © . 

Another assumption used in our symmetry consider- 
ations (Section [TTJ) and elsewhere throughout the paper 
is that the diffracted wave is plane. In fact, the wave 
diffracted by a small particle is spherical at a large dis- 
tance from the latter. But on a scale much smaller 
than the distance between the particle and the measuring 
point the wave may be considered planed. Its intensity, 
indeed, depends on the distance from the particle. 

Since the experiments on light diffraction by the lattice 
of magnetic vortices have been carried out for the linear 
polarizations of incident light, we also use these polar- 
izations in theoretical calculations. The following desig- 
nations for the linear polarizations are used throughout 
the paper: the wave polarized so that its electric field 
vector lies in the plane defined by the k vector and the 
axis of the magnetic vortex e z is termed p-polarized (see 
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Fig. [IJ; the wave with the electric field perpendicular to 
that plane is s-polarized. 

One more important approximation consists in re- 
stricting ourselves to the first order in the magnetic mo- 
ment, which mathematically corresponds to the first or- 
der in a small parameter j/e (see Eq. (J7J). Our choice 
of approximation is determined by both the fact that 
the phenomenon described by Eq. (jj) is linear in M 
and that 7/e is a very small parameter for the existing 
ferromagnet o 32 ! 33 . 



IV. THE BORN APPROXIMATION 

The simplest calculation is based on the assumption 
that permittivity tensor e is almost equal to unit tensor 
1: 



e = 1 + Se. 



(9) 



The Maxwell equations can be solved in this case for a 
scattering particle of an arbitrary shape 31 . According to 
the well-known theory, the electric field of diffracted wave 
E SC at can be written as 



E s 



ik'Ro 



4ttR 



k' X 


k'x / 




iv 



<kE e lqr dV 



(10) 



Here Eo is the electric field of the incident wave, q stands 
for k — k', V is the volume of a scattering particle, Rq is 
the distance between the scattering particle and the point 
where E sca t is measured. According to our assumptions 
f Section ing . Rq is much bigger than the linear dimension 
of scatterer V 1 ^ 3 . 

We can now use ([7]) and © to determine Se and cal- 
culate the diffracted electric field. It is convenient to use 
the matrix form to represent the diffraction coefficients 
for different linear polarizations: 



g _ ^0 f S ss S, 



sp 



(11) 



here the diagonal terms represent the diffraction with- 
out a change of polarization, while the off-diagonal terms 
stand for the diffraction to another polarization. The 
term k$ in Eq. (fTTI) is the absolute value of the wave 
vector fco — uj/c. The matrix components in are 



S, 



(e - 1) e lq W 



S sp = i-f(k- J M (r)e iqr dV 



(12) 



(13) 



S ps = -i 7 [k'- y^MWe^dV) (14) 



S pp = (k • k') f (e - 1) e iq W - 
Jv 

[k x k'] ■ / M (r)e iqr dV 



-fry 



(15) 



When analyzing Eqs (fT^ - (Tl"5l) we should first of all 
note that the linear polarizations do not mix without 
magnetization. Since we are interested in the inten- 
sity effect that is linear in M, only the diagonal terms 
of matrix (JTTJ) should be taken into account. The ma- 
trix component S ss does not depend on the magnetic 
moment, hence, we come to a conclusion that the s- 
polarized light diffraction does not feature non-reciprocal 
properties, while for the p-polarized light there is a non- 
reciprocal contribution in the diffraction intensity. 

Finally, taking into account Eq. (|5|), we can expand 
the term e lqr into Fourier series and neglect all but the 
first two terms (e lqr * 1 + iqr). These correspond to 
the electrodipole and electroquadrupole terms that are 
of the zero and the first order in a/ A. Calculation of the 
intensity for M as given by ([6]) leads to 



I, = 



ml 



1 



(16) 



" 9iZg 
3koa 
"LIT 



■(|e-l| 2 (n-n') 2 -((n + n')-e z )x 
((n-n')-(n.n') 2 )i?e(( £ -l)* 7 )) (17) 



Here I p and I s are the intensities of the diffracted light 
for p and s polarizations of incident light, n and n' are 
the wave vectors of the incident and diffracted waves nor- 
malized to unity: n = k/ko, n' = k'/ko, e z is the unit 
vector along the vortex axis. 

In order to find the intensity of diffracted light for 
the unpolarized incident light (we label it I np ), the in- 
tensity should be averaged over the polarizations. For- 
mally, it leads to I np = \ (I s + I p ), where I np ~ a. From 
Eqs (Tl6l fT7]) it is clear that the above calculation is con- 
sistent with the phenomenological considerations. 

In the conclusion of this section, we point out that 
the non-reciprocal light diffraction was of the same or- 
der of magnitude for both s and p polarizations of in- 
cident light in the experiment 2 ^. But within the used 
approach the non-reciprocal contribution exists only for 
the p-polarized light. In the next section we propose a 
more general theory that has no such shortcoming. 



V. PERTURBATION THEORY 

A more exact solution of the problem is gained through 
the perturbation theory based on condition (|5|). To 
account for the electroquadrupole and magnetodipolc 
terms, we need to expand the solution up to the first 
order in a/ A. The correction to the electrodipole term 
will be gained naturally. Unlike in some simple models^ 
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we cannot use the quasistatic Maxwell equations here. 
Hence, our approach is based on the approximations of 
the Maxwell equations 



rotE = ik H, 
rotH = — ifc eE. 
div (eE) = 0, 
divH = 0, 

that follow from the estimation 

rotE ~ E/a, 
rotU ~ H/a. 



(18) 
(19) 
(20) 
(21) 

(22) 
(23) 



Next, we assume that the permittivity tensor is not large. 
Taking into account that j/e << 1, we impose the condi- 
tion one: |e| << A/a (A/a >> 1 due to ([5])). Introducing 
a typical scale of fields in the particle L ~ A/-%/|e| (for 
real e L is the wavelength in the medium, for imaginary 
e it is the skin depth), we have 



a « L 2 /X. 



(24) 



This condition along with Eq. (O allows us to derive the 
equations for the fields of zero and first order in a/ A from 
Eqs (HEM]). 

First of all, let us find the electric field inside the par- 
ticle. As is well-known, in zero order in a/ A and in "f/e 
the electric and magnetic fields have the form 



0,0 



:Eo, 



E' 

£ + 2 

H°<° = Ho = [n x E ] 



(25) 
(26) 



The first upper index in E 00 and H°'° is used to denote 
the order in a/A, the second index specifies the order in 
7/e. 

The equation for E 1 - follows from Eqs (0 and (l2"2l 

rotE 1 ' = ik H°'°. (27) 

By solving this equation with a zero boundary condition 
at the infinity we arrive at 



E 1 ' = -((k-r)Eo-(E -r)k) + 
((k-r)E + (E -r)k). 



2 2e + 3 



(28) 



We now calculate the fields linear in the magnetic mo- 
ment. The equation for H 0,1 is 



rotH. ' 1 = 0, 
divH - 1 = 0, 



(29) 
(30) 



and, taking into account the zero boundary condition at 
the infinity, we have H° 1 = 0. 

Finally, the equations for E 01 and E 1,1 are very sim- 
ilar: 

roiE}' 1 = 0, (31) 
dro (eE 1 - 1 + i 7 [E 1 ' xM(r)]) = 0, (32) 



where I = 0, 1. Using Eq. (|3"Tj) . the solution can be found 
as E 1,1 = — Vip 1 (r), where ip is the scalar function that 
satisfies the equation 



At/) 1 = i-div [E 1 - x M] 



(33) 



The right-hand part of Eq. (133)) may be treated as a 
charge density distribution. Solving this equation we 
have 



i 7 

47T £ 



div [E 1 - (r') x M (r')] 



dV. 



(34) 



Substitution of Eqs $25$, dHJ), @ in Eq. yields a 
formal solution to the problem. Nevertheless, it should be 
mentioned that the integral over the particle in the right- 
hand part of Eq. (|3"4")l cannot be expressed as elementary 
functions. 

The next step of our solution procedure is calculation 
of the electric current density j that is related to the 
electric field^: 



J = -4^( e - 1 ) K 
Substitution of ([7]) in (f3"5"|) yields 



:i,o 



ICJ 



r- = -^(e-DE 



1.0 



(35) 



(36) 



jU = -™(e- i) e 1 ' 1 + ^— \E X - Q x Ml . (37) 

47T 47T 

Knowing the electric current density distribution over 
the particle, we can calculate the electric field of the 
diffracted wave far from the diffracting particle^: 

E scat = -2-— / j (r) e- ik r dV. (38) 

We expand the expression e _,k r in this formula in the 
Fourier series and keep the first two terms (e _lk r ~ 
1 — ik'r). This expansion corresponds to the inclusion 
of the electroquadrupole, magnetodipole and the first- 
order correction to the electrodipole term. All of these 
terms have the same order of magnitude in a/ A. The 
waves radiated by them all interfere with the zero-order 
electrodipole radiation. There is no point in separating 
them. 

Finally, for the electric field of the diffracted wave we 
have 



E 



scat " 
0.1 



JJ^gikoRo-iwt 



(0°'° 



■j 1 ' 1 )-i(r.k')(i ' +j " 1 ))dV; (39) 



where j 0,0 , j 1,0 , j 0,1 , j 1,1 are determined by Eqs (f3"61 |3"T1) . 
We have earlier noted that the integral defining the elec- 
tric field inside the particle cannot be expressed as el- 
ementary functions; however, the double integral that 
appears in Eq. (|39l) is taken by changing the integration 
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order. Thus, E sca t turns out to be expressed as elemen- 
tary functions, indeed. 

The last step of the calculation is to find the intensity 
of the diffracted light. Again, for two linear polarizations 
of the incident light we get 



I, = 



1 



e + 2 



(l-((n 
irkoa 



n 



32 



Re 



(r, 



2e + 3 



(40) 



fc 4 q 6 



e-1 



e + 2 
irkoa 



(n-n') -((n + n')-e z )x 



32 



67 



(n ■ n') + 



2e + 3 
x nn • n') — (n • n' 



s- 1 

\2' 



.(41) 



Eqs (|40l 141]) manifest that the non-reciprocal contribu- 
tion exists for both s- and p-polarized incident light. 
For the p-polarized light they show a somewhat differ- 
ent dependence on the angle between n and n' from that 
yielded by the Born approximation. These formulas are 
completely consistent with the phenomenological consid- 
erations f Section HIl). 

The perturbation theory that accurately takes into ac- 
count all of the terms linear in a/ A qualitatively agrees 
with the experiment 2 -. To compare the results of our the- 
oretical calculation with the experimental data we per- 
formed simple estimations for the parameters of cobalt^. 
The wavelength is A = 632. 8nm. At this wavelength 
cobalt absorbs light quite well, so the permittivity e is 
defined by the complex constants e — — 12.6 + 22.88i, 7 = 
0.749 — 0.602i. In the experiment 26 the particles are fiat 
(their thickness is much less than their lateral size). The 
scale of variation of the internal current density that 
is important for our calculation is a particle thickness. 
So we take the particle radius a = 30nm. The inci- 
dent light propagates along the vortex axis (n = e z ) 
and the diffracted light is deflected from it by 30°. For 



these parameters we have AI S /I S « 1.12- 10 -2 , AI p /I p 
1.29-10 -2 . The experimental value of the effect is 2-10~ 3 . 
So, we come to a conclusion that although this estima- 
tion is very rough it turns out to be a satisfactory fit to 
the experiment. The difference is apparently attributed 
to the different shape of the particle and to the violation 
of the assumption that the permittivity value is not large 
(see (ED- 



VI. CONCLUSION 

We have solved the problem of light diffraction by a 
spheric magnetic particle with the vortex magnetization 
distribution up to the first order in the a/ A parameter. 
It should be noted that the approach can easily be ex- 
panded to the next-order calculation, but generally it is 
not necessary in order to gain an insight into the na- 
ture of these phenomena. Our results show that the 
vorticity-dcpcndcnt non-reciprocal contribution in the in- 
tensity of the diffracted light takes place for both s and 
p polarizations of the incident light. The non-reciprocal 
term arises due to the interference of the zero-order elec- 
trodipole radiation and the first-order electrodipole, elec- 
troquadrupole and magnetic dipole radiation that lin- 
early depends on the magnetic moment of the particle. 
The estimations of the effect value for the parameters of 
cobalt yield a relative value AI/I ~ 10~ 2 that fits the 
experiment^ by an order of magnitude. 
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